function intg = numint(So,R,T)